Mon. Not. R. Astron. Soc. 000, 1-?? (2007) Printed 2 February 2008 (MN WF^ style file v2.2) 



QSO-LRG 2-Point Cross-Correlation Function and 
Redshift- Space Distortions 

. G. Mountrichas/'*', T. Shanks^, S. M. Croom^, U. Sawangwit^, D. P. Schneider^, 
A. D. Myers^ K. Pimbblet^ 

^Department of Physics, University of Durham, South Road, Durham DHl 3LE, UK 
^ School of Physics, University of Sydney, NSW 2006 Australia 

^Department of Astronomy and Astrophyscics, Penn State, 504 Davey Laboratory, University Park, Pennsylvania 16802 
, '^Department of Astronomy, MC-221, University of Illinois, 1002 W. Green Street, Urbana, IL 61801 USA 
■ ^Department of Physics, University of Queensland, Brisbane, Qld 4072, Australia 



27 March 2007 



ABSTRACT 

We have measured the bias of QSOs as a function of QSO luminosity at fixed 
redshift {z < 1) by cross-correlating them with Luminous Red Galaxies (LRGs) in 
the same spatial volume, hence breaking the degeneracy between QSO luminosity 
and redshift. We use three QSO samples from 2SLAQ, 2QZ and SDSS covering a 
QSO absolute magnitude range, —24.5 < Alf,, < —21.5, and cross-correlate them with 
2SLAQ {z « 0.5) and AAOmega {z « 0.7) photometric and spectroscopic LRGs in the 
same redshift ranges. The spectroscopic QSO samples, in the spectroscopic LRG areas, 
contain 300-700 QSOs and in photometric LRG areas up to 7000 QSOs. The 2-D and 3- 
D cross-clustering measurements are generally in good agreement. Our (2SLAQ) QSO- 
LRG clustering amplitude (rg — G.Slg'gh^^Mpc) as measured from the semi-projected 
cross-correlation function appears similar to the (2SLAQ) LRG-LRG auto-correlation 
amplitude (ro = 7.45 ± 0.35h-iMpc) and both are higher than the (2QZ-h2SLAQ) 
QSO-QSO amplitude (ro ~ 5.0h~^Mpc). Our measurements show remarkably little 
QSO-LRG cross-clustering dependence on QSO luminosity. If anything, the results 
imply that brighter QSOs may be less highly biased than faint QSOs, the opposite 
direction expected from simple high peaks biasing models, where more luminous QSOs 
are assumed to occupy rarer high mass peaks. Assuming a standard ACDM model and 
values for bLuc measured from LRG autocorrelation analyses, we find bq = 1.45 ±0.11 
at Mbj « -24 and bq = 1.90 ± 0.16 at Affc, « -22. 

We also find consistent results for the QSO bias from a z— space distortion analysis 
of the QSO-LRG cross-clustering at z « 0.55. The velocity dispersions fitted to QSO- 
LRG cross-correlation, ^(cr, tt), at ±680 kms~^ are intermediate between those for 
QSO-QSO and LRG-LRG clustering, as expected given the larger QSO redshift errors. 
The dynamical infall results give /3q = 0.55±0.10, implying bg = 1.4±0.2. Thus both 
the z— space distortion and the amplitude analyses yield 6q « 1.5 at M^j « —23. The 
implied dark matter halo mass inhabited by QSOs at z « 0.55 is ^ 10^^/i~^Mq, again 
approximately independent of QSO luminosity. 

Key words: galaxies:clusters:general-quasars:general-cosmology:observations-large- 
scale structure of Universe 



1 INTRODUCTION 

Useful constraints on the nature of QSOs can be drawn from 
even the simplest measure of QSO clustering, the ampli- 
tude of the real-space 2-point correlation function. For ex- 
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ample, measuring QSO clustering on small scales can con- 
strain models of galaxy mergers and quasar formation (My- 
ers et al. 2007). Furthermore, the redshift evolution and lu- 
minosity dependence of the QSO bias can be studied. Recent 
work on the clustering dependence of the QSO bias on red- 
shift suggested evolution of the 2QZ QSO bias (Groom et 
al. 2005). However, the fact that the most luminous QSOs 
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lie at high redshifls, while the faintest are at low redshifts, 
made it difficult to study how QSO properties depend on lu- 
minosity. Now, surveys such as 2SLAQ (Cannon et al. 2006) 
that span a wide range of QSO luminosities, have broken 
that rcdshift-luminosity degeneracy and revealed little QSO 
clustering dependence on luminosity, at fixed redshift (da 
yingela et al. 2006). Moreover, the amplitude of the QSO 
clustering is correlated with the average mass of the halos 
associated with the QSOs, providing indications of QSO life- 
times and making it possible to constrain QSO evolutionary 
models (Groom et al. 2005, da Angela et al. 2008). 

Redshift space distortions of the clustering pattern also 
contain dynamical information on the QSO bias that arc 
independent of any assumption about underlying mass clus- 
tering. The clustering is affected at small scales by the rms 
velocity dispersion of QSOs along the line-of-sight (Fingers 
of god) and by dynamical infall of matter into higher density 
regions, which causes a flattening of the clustering pattern in 
the redshift direction. In addition to these dynamical distor- 
tions, geometrical distortions are introduced if an incorrect 
cosmological model is used in order to convert the observed 
redshifts into comoving distances. They therefore also con- 
strain, more weakly, the value of the cosmological density 
parameter, il,rn- 

In the linear regime of clustering, dynamical infall is 
governed by the parameter /3 = f2^®/&. (Peebles, 1980, 
Kaiser 1987, Loveday et al. 1996, Matsubara & Suto 1996, 
Matsubara & Szalay 2001, Peacock ct al. 2001, Hoylo ct 
al. 2002, Coil et al. 2005, Myers et al. 2006, 2007, Por- 
ciani & Norberg 2006, da Angela et al. 2008, Ross et al. 
2007). In recent years, measurements of QSO clustering (da 
Angela et al. 2008) yielded a i3qso{z = 1.4) = O.GOlJ];" 
and bQso{z = 1.4) — 1.5 ± 0.2 for a combined sample of 
QSOs from the 2dF QSO Redshift Survey (2QZ, Croom et 
al. 2004) and the 2dF-SDSS LRG and QSO Survey (2SLAQ, 
Cannon et al. 2006). Ross et al. (2007) performed simi- 
lar measurements on the 2SLAQ Luminous Red Galaxies 
(LRGs) clustering and found I3lrg{z = 0.55) = 0.45lo:o5 
and bLRoiz = 0.55) = 1.66 ± 0.35. 

In this paper we use QSOs from the 2QZ, 2SLAQ and 
SDSS (York ct al. 2000) Data Release 5 (DR5; Adelman- 
McCarthy 2007) surveys and LRGs from 2SLAQ and 
AAOmega first to study the dependence of QSO-LRG clus- 
tering amplitude on QSO luminosity. We also measure QSO- 
LRG redshift distortions to estimate the dynamical infall pa- 
rameter, /3. These surveys provide largo numbers of QSOs 
and LRGs with a range of luminosities at fixed redshifts. So 
our results, for example on QSO bias, should be statistically 
improved over those from QSO-QSO clustering. Moreover, 
we use photometric LRG samples, which are significantly 
larger than the spectroscopic ones and measure QSO 2-D 
correlation function amplitudes, using Limber's formula to 
convert to 3-D real-space measurements. 

In section 2 we describe the QSO and LRG samples we 
use in our measurements. In section 3 we present results from 
the 2-point cross-correlation function cij{6), in Sections 4 and 
5 we show our results for the redshift-space cross-correlation 
function, ^s, and the semi-projected cross-correlation func- 
tion, Wp(a)/a, respectively. Section 6 has our results for 
the real-space cross-correlation function and in Section 7 we 
present our ^(cr, tt) results, model the redshift-space distor- 
tions and estimate the QSO infall mass and bias. We also 



Table 1. The numbers of spectroscopic QSOs and spectroscopic 
2SLAQ LRGs. 



2SLAQ area (0.35 ^ ^ ^ 0.75) 




QSOs 


LRGs 


2SLAQ 


699 


8,656 


2QZ 


307 


5,995 


SDSS 


218 


5,995 



Table 2. The numbers of spectroscopic QSOs and photometric 
LRGs. 





2SLAQ 


area (0.35 ^ z 


< 0.75) AAOmega 


area (0.45 < z < 0.90) 




QSOs 


LRGs 


QSOs 


LRGs 


2SLAQ 


503 


19, 300 


786 


23, 836 


2QZ 


1,048 


32, 188 


1,265 


40, 060 


SDSS 


7,395 


468, 416 


7,083 


571, 676 



study, in Section 8, the dependence of the redshift-space 
cross-correlation function, the QSO bias and the mass of 
the QSO-LRG Dark Matter Haloes (Mdmh) on QSO lumi- 
nosities at fixed redshifts. Finally, in Section 9 we present 
our conclusions. 



2 DATA 

2.1 Spectroscopic data 

Our spectroscopic LRG sample is taken from the 2SLAQ 
and consists of 8,656 LRGs within 0.35 ^ z ^ 0.75 (that 
is the 'Gold Sample' of Ross et al. (2007), see their Fig. 1 
for the LRG distribution). 5,995 LRGs are in the Northern 
strip (sectors a, b, c, d, e, from Fig. 1 of Ross et al. 2007) 
and 2,661 LRGs in the Southern strip (sector s). 

Our QSOs are taken from three spectroscopic samples; 
the NGC of 2QZ, the NGC-FSGC of 2SLAQ and the Data 
Release 5 (DR5) of SDSS. The QSO redshift range mainly 
used in our analysis is the same one as for the 2SLAQ LRGs 
(0.35 < 2 «; 0.75). The 2SLAQ QSOs have the same distri- 
bution on the sky as the spectroscopic 2SLAQ LRGs (see 
Fig. 2 of da ingela et al. 2008). The 2QZ and SDSS QSOs 
cover only the NGC of the 2SLAQ LRGs. The brightest of 
our QSO samples is from SDSS, which consists of QSOs 
with iAB < 19.1. Our 2QZ sample consists of QSOs with 
18.25 ^ 6j ^ 20.85 and the 2SLAQ sample of QSOs with 
18.0 ^ gr < 21.85. After matching the QSO and spectro- 
scopic 2SLAQ LRG areas we get the numbers shown in Ta- 
ble 1. 

2.2 Photometric data 

For measuring the 2-point angular cross-correlation func- 
tion, w(9), we can also use larger, photometric LRG sam- 
ples from the 2SLAQ and the AAOmega surveys. For the 
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AAOmega LRGs wc have followed the selection criteria of 
Ross ct al. (2007). The numbers for these LRG data sets 
and the new matched QSO samples axe shown in Table 2. 
We should mention here that in the case of photometric 
2SLAQ and AAOmega LRGs with the 2SLAQ QSOs we 
did not perfectly match the two areas. The photometric 
2SLAQ LRG area is a rectangle (137.5° ^ ra ^ 230.0°, 
-1.25° ^ dec < 1.0°) which covers the whole NGC of 
2SLAQ, including the gaps in between the sectors a, b, c, 
d, o, whereas the 2SLAQ QSOs have the distribution men- 
tioned in the previous Section. This slight mis-match does 
not affect our results, but is the reason that, in Table 2, 
photometric 2SLAQ LRGs matched with 2SLAQ QSOs ap- 
pear more rmmerous than the spectroscopic 2SLAQ LRGs 
matched with 2SLAQ QSOs in Table 1 (19,300 vs. 8,656). 
Furthermore, the 2SLAQ QSO set, in the 2SLAQ photomet- 
ric area, is smaller than that in the 2SLAQ spectroscopic 
area (503 vs 699) because in the photometric case we do not 
use the SGC part of 2SLAQ (sector s). 

The redshift distribution of all the QSO samples, in a 
redshift range of 0.0 < a < 1.0 is shown in Fig. 1. The 
2SLAQ and 2QZ distributions are, as expected, very similar. 
The SDSS distribution is roughly flat, in the redshift range 
we are interested in, although we note a peak at z ~ 0.2 
which is outside our region of interest. 



3 QSO-LRG ANGULAR 

CROSS-CORRELATION FUNCTION 

3.1 Cross-correlation and error estimators 

In this section we first estimate the 2SLAQ QSO— 2SLAQ 
LRG and 2QZ QSO-2SLAQ LRG angular cross-correlation 
functions w(8) using our spectroscopic samples, both for 
QSOs and LRGs. We then use our photometric 2SLAQ 
and AAOmega LRG samples and cross-correlate them with 
2SLAQ, 2QZ and SDSS (DR5) QSOs. The estimator we use, 
in all cases, is 



w{e) 



DD{9) Nrd 
DR{e) Nlrg 



(1) 



where Nrd is the number of points in our random catalogue, 
Nlrg is the number of LRGs, DD(0) are the data-data 
pairs, i.e. QSO-LRG pairs and DR{6) are the QSO-random 
point pairs counted at angular separation 0. 

Our spectroscopic LRG random catalogue is the same as 
that described by Ross et al. (2007) (with 20 x more randoms 
than LRGs) whereas the photometric 2SLAQ and AAOmega 
LRG random catalogues contain ~ llx more randoms than 
LRGs. Coverage completeness and spectroscopic complete- 
ness are taken into account in the construction of all random 
catalogues. 

The error estimator we use throughout this paper is the 
so-called Field-to-Field estimator. The accuracy of the errors 

on our measurements plays an important role, particularly 
when we use our results to model the redshift-space distor- 
tions in Section 8. The reason for not using a different error 

estimator, such as the Poisson estimator (o'(6) = -^^^p^) 
is that it becomes increasingly inaccurate at larger scales, 
as the QSO-LRG pairs become less independent. In order 
to calculate the Field-to-Field errors we have divided the 



2SLAQ area into 4 approximately equal areas and then mea- 
sure the cross-correlation functions in each one of these ar- 
eas. The Field-to-Field error is then given by the following 
expression (Myers et al. 2005) 



N 



(2) 



where DRl{9) is the data-random pairs in the subarea, 
DR{9) is the overall number of data-random pairs, u>L{d) is 
the correlation function measured in the subarea and u){9) 
is the overall correlation function. 



3.2 w{0) results from the redshift samples and 
correction for fibre collision 

Fig. 2 shows our angular correlation results when wo cross- 
correlate 2SLAQ QSOs with spectroscopic 2SLAQ LRGs. 
Filled circles show the results when we use our whole 2SLAQ 
QSO sample and cross-correlate it with our sample 8 2SLAQ 
LRGs. Open circles show the results from 2SLAQ QSOs 
in the redshift range of 1.0 < z < 2.2 cross-correlated 
with the same LRG sample. Finally, triangles show the re- 
sults when we use 2SLAQ QSOs in the redshift range of 
0.35 ^ z ^ 0.75. The fibre collision problem affects our re- 
sults as can be seen from the anti-correlation we find, in all 
cases, between our QSOs and LRGs on small scales. The 
results are very similar when we use 2SLAQ QSOs in the 
whole redshift range and within 1.0 ^ « ^ 2.2. The anti- 
correlation amplitude is increased when we use our low red- 
shift QSO sample, which also shows an increased fibre ef- 
fect, but its statistical significance is smaller as the sample 
is much smaller than the other two. Now, we shall use these 
results in order to correct our 2-D and 3-D results, for fibre 
collisions. 

We shall base our correction on the Hawkins et al. 
(2003) expression, i.e. 

1 + Wp 

™/ = T— {i) 

1 + 

Wp is the cross-correlation results when our samples come 
from the input catalogue; in our QSO-LRG case = if 
there is no redshift overlap between QSOs and LRGs and 
assuming no lensing. represents the cross-correlation re- 
sults from the 2SLAQ catalogue and those arc the results 
shown in Fig. 2. The correction should therefore be 



Wf = 



(4) 



1 + w^ 

We determine Wz based on the measurements shown 
in Fig. 2. As already noted, on small scales (^ 2'), we 
see the anti-correlation caused by the fibre effect. We also 
note a small positive bump on scales up to ~ 80'. This 
bump is expected to be caused by physical QSO-LRG cross- 
correlations for the 0.35 ^ z ^ 0.75 range (triangles) and 
also to a lesser extent in the unrestricted redshift range 
(closed circles). The cause of small positive excess seen at 
the 1.0 < z < 2.2 could be due to the fibre collision effect. 
If we include this small positive bump in the fibre correc- 
tion, our ^(s) cross-correlation measurements do not change; 
therefore this bump will be ignored when we present our 3-D 
results in the next Section. Nevertheless, the feature's effect 
is larger for the 2SLAQ QSO-2SLAQ (photometric) LRG 
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Figure 1. The solid line shows N(z) for the 2SLAQ sample, the dashed line for the 2QZ and the dotted line for the SDSS sample. The 
long dashed line and the dashed-dotted line show the redshift range we use for our measurements, when we cross-correlate our QSO 
samples with the 2SLAQ and AAOmega LRGs, respectively. The 2SLAQ and 2QZ distributions are, as expected, very similar. The SDSS 
distribution is flat, in the 2SLAQ and AAOmega redshift ranges we are interested in. 



w{6) measurements which will be presented later (Fig. 4) 
and it will be taken into account there. 

Finally, Fig. 3 shows our results when wo cross-correlate 
2QZ QSOs with spectroscopic 2SLAQ LRGs. No fibre col- 
lision effect is expected in this case since 2QZ QSOs had 
higher priority than 2dFGRS galaxies for spectroscopic ob- 
servations. Thus there is no issue of fibre incompleteness 
(see also, e.g., Myers et al. 2003). The reason for perform- 
ing these measurements is to see if we can detect any signal 
in the case of the overlapped redshift range (triangles) . Al- 
though the results when we use 2QZ QSOs in 1.0 ^ z ^ 2.2 
range (open circles) seem to give a null average signal as ex- 
pected, in the case of the common QSO-LRG redshift range 
(triangles) the samples are small and the results appear too 
noisy to draw any statistically significant conclusions. 



3.3 Results from the photometric samples 

Before we proceed to measure the rcdshift-space cross- 
correlation function, we shall repeat our w{9) measurements, 
cross-correlating our spectroscopic QSO samples from SDSS, 
2QZ and 2SLAQ with the photometric LRG samples, which 
were presented in Section 2.2 and Table 2. The purpose of 
these measurements is to see how the LRGs correlate with 
bright and faint QSOs and then use Limber's formula to con- 
vert these 2-D measurements into 3-D real-space measure- 
ments and compare them with our results from the spectro- 
scopic samples (Section 4). 

The results using 2SLAQ photometric LRGs are shown 
in Fig. 4 (filled circles axe corrected for fibre collisions, from 
0.1' < 61 < 100', i.e. including the bump in Fig. 2) and the 



results using the AAOmega photometric LRGs are shown in 
Fig. 5. We should note that neither the AAOmega nor the 
2SLAQ results have been corrected for stellar contamination 
of the LRGs (~ 15% in AAOmega and ^ 5% in 2SLAQ). 
In both cases, triangles show the results using SDSS QSOs, 
filled circles using 2SLAQ QSOs and open circles using 2QZ 
QSOs. The 9{] and ro values from the fits (in the range 0.1'- 
100' ) to these measurements are also shown. Both results 
show a slightly steeper slope for the brightest QSO sample 
(SDSS, 7 = -0.8±0.1) compared to the other QSO samples. 

Comparing our w{6) results from 2QZ QSO- 
spectroscopic LRGs (triangles of Fig. 3) with those 
from 2QZ QSO-photometric LRGs (open circles of Fig. 4) 
we note that the latter shows a positive signal whereas the 
former shows no signal, at least on small scales, and appears 
to be noisy. This is probably because the 2QZ QSO sample 
is much larger (« 3 x ) in the photometric LRG area (Table 
2) than it is in the spectroscopic LRG area (Table 1). 



4 3-D CROSS-CORRELATION FUNCTIONS, 

as) AND e(r) 

In this Section we present our results for the QSO-LRG 
redshift-space cross- correlation function, ^(s), for different 
QSO and LRG samples. Fig. 6 shows the results using 
our 2SLAQ QS0-2SLAQ (spectroscopic) LRG samples. Our 
2SLAQ QSO sample consists of QSOs with the same aver- 
age redshift as the 2SLAQ LRGs, i.e. 0.35 < ^ ^ 0.75. Filled 
circles show the results when the fibre collision effect is not 
taken into account and open circles show the results when 
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Figure 2. 2SLAQ QSO-2SLAQ (spectroscopic) LRG cross-correlation results. Filled circles show the results when we use our entire 
(9,044) 2SLAQ QSO sample and cross-correlate it with our spectroscopic 2SLAQ LRGs (9,856). Open circles show the results when we 
use (6,002) 2SLAQ QSOs in the redshift range of 1.0 ^ x ^ 2.2 and finally, triangles show the results when we use (699) 2SLAQ QSOs 
with a redshift range of 0.35 ^ z ^ 0.75. The fibre collision effect is seen in the results as we get an anti-correlation between our QSOs 
and LRGs. The solid line shows our best fit to the results (filled circles). 
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Figure 3. 2QZ QSO-2SLAQ LRG cross-correlation results. Filled circles show the results when we use our entire (2,854) 2QZ QSO 
sample ('11' quality, e.g.. Groom et al. 2004) and cross-correlate it with our spectroscopic 2SLAQ LRGs (NGG, 5,995 LRGs). Open 
circles show the results when we use (1,699) 2QZ QSOs in the redshift range of 1.0 ^ z ^ 2.2 cross-correlated with the same LRG 
sample. Finally, triangles show the results when we use (307) 2QZ QSOs with redshift range 0.35 ^ z ^ 0.75. 
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Figure 4. Triangles show the results from SDSS QSOs, filled circles from 2SLAQ QSOs (corrected for fibre collisions) and open circles 
from 2QZ QSOs cross-correlated with photometric 2SLAQ LRGs. We have also plotted the fits to these measurements. We note that the 
slope is steeper for the brightest SDSS QSO sample comparing to the fainter 2QZ and 2SLAQ QSO samples. 
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Figure 5. Triangles show the results from SDSS QSOs, filled circles from 2SLAQ QSOs and open circles from 2QZ QSOs cross-correlated 
with photometric AAOmega LRGs. We have also plotted the fits to these measurements. Once more, we note that the slope is steeper 
for the brightest SDSS QSO sample comparing to the fainter 2QZ and 2SLAQ QSO samples. 
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the fibre collision is taken into consideration, as described 
in the previous Section. As wo can sec the effect is bigger on 
small scales (i.e. s < 3h~^Mpc) than it is on larger scales. 

We also repeat this redshift-space cross-correlation mea- 
surement using the 2QZ and the SDSS QSO samples, both 
within a redshift range of 0.35 ^ « < 0.75, and spectroscopic 
2SLAQ LRGs. As already mentioned, the 2QZ survey has 
a brighter magnitude limit than 2SLAQ, i.e. bj = 20.85 
instead of bj = 21.85. The brightest of the three samples 
is that from SDSS (jab < 19.1). Our results are shown in 
Fig. 7. On small scales, ^ 5h~^Mpc we can see the suppres- 
sion of f (s), due to the non-linear redshift-space distortions 
(small-scale peculiar velocities, {w'^)^^^). The effect is loss 
for SDSS, possibly because SDSS uses narrow linos to de- 
termine rcdshifts for QSOs at low redshift and [OIII] 5007 
is seen for the whole redshift range in question for SDSS 
spectra. 2SLAQ/2QZ redshifts are purely template fits and 
so will instead bo driven by the broad lines. On larger scales 
the results arc in good statistical agreement, regardless of 
the brightness of the QSO sample. This agreement is also 
confirmed by the fits to these measurements (Table 4). Since 
the results are affected by the Finger of god effect on small 
scales, the fits arc applied on 5 — 25h^^Mpc scales. The 
agreement between the ^(s) results, suggests that QSO bias 
is also independent of QSO luminosity since the QSOs span 
a range of more than 2 magnitudes at fixed redshift. 

To further check our observations, we use Limber's for- 
mula (Limber 1953, Rubin 1954) and following Phillipps ct 
al. (1978) we calculate the real-space cross-correlation func- 
tion, ^(r) from our previous (Section 3) w{9) measurements 
from the photometric LRG samples. The fits appear in Fig. 
8 and in Table 3. The black solid line (ro = 7.5 ± 0.3, 7 = 
1.7 ± 0.2) is the fit using the 2SLAQ QS0-2SLAQ LRG 
sample, the dotted line (ro = 8.0 ± 0.4, 7 = 1.7 ± 0.2) us- 
ing the 2QZ QS0-2SLAQ LRG sample and the dashed Une 
(ro = 7.0 ± 0.3, 7 = 1.8 ± 0.1) using the SDSS QSO-2SLAQ 
LRG sample. The blue lines are the fits using our QSO sam- 
ples with AAOmega LRGs. Once again, we note that the 
results are approximately independent of the luminosity of 
the QSO sample. 

In Fig. 7 we have plotted the fits from the QS0-2SLAQ 
(photometric) LRGs with the results from the spectroscopic 
2SLAQ LRGs already discussed. This is a comparison be- 
tween ^(r) (photometric) and ^(s) (spectroscopic) results. 
As we can see from Fig. 7 the agreement is not good at small 
scales (< 5h^^Mpc) but this is duo to non-linear redshift- 
space distortions (Finger of god) that affect ^(s) but not 
the ^(r) measurements on small scales. A fairer comparison 
can be made on larger scales. Taking into account the Kaiser 
boost (~ 1.25 for /? = 0.35, as we shall sec next) we compare 
the fits (Tables 3 and 4) to the QS0-2SLAQ (photometric) 
LRGs and QSO-2SLAQ (spectroscopic) LRGs, respectively. 
We see that they are in reasonable agreement. 

Finally, in Fig. 9 we compare our redshift-space mea- 
surements for QSO-LRGs with QSO-QSO (da Angela et al. 
2008) and LRG-LRG (Ross et al. 2007) measurements. Tri- 
angles show the 2QZ-I-2SLAQ QSO {{s) results, filled cir- 
cles show the 2SLAQ LRG ^(s) results and open circles 
show our results for the 2SLAQ QS0-2SLAQ (spectroscopic) 
LRG redshift-space cross-correlation. The fact that QSO- 
QSO and QSO-LRG results appear flatter than the LRG- 
LRG ones, on small scales, may be explained by the in- 



trinsic dispersion and the redshift errors, which are higher 
for broad emission line QSOs than for LRGs (~ 650kms^^ 
vs. w 300kms"^). This would affect mostly the QSO-QSO 
results and the least the LRG-LRG results, as observed. 
The QSO-LRG measurements should then lie in between the 
QSO-QSO and LRG-LRG measurements (as they do). On 
larger scales (^ 5h~^Mpc) the QSO-QSO correlation am- 
plitude appears lower than the LRG-LRG one. This implies 
that the QSO bias is smaller than the LRG bias. Assum- 
ing the model ^mm = -^^^ (see Section 7.2), the agreement 
between the QSO-LRG amplitude and the LRG-LRG am- 
plitude would imply that bg ^ bL- So the overall conclusion 
is that 6q < 6i. These issues will be further discussed in 
Section 8. 



5 THE SEMI-PROJECTED 

CROSS-CORRELATION FUNCTION 

If si and S2 are the distances of two objects 1, 2, measured in 
redshift-space, and 6 the angular separation between them, 

then a and tt are defined as 

TT = (s2 — si), along the line-of-sight (5) 

a = -^—^—^9, across the line-of-sight (6) 

The effects of redshift distortion appear only in the radial 
component, tt; integrating along the tt direction, we calculate 
what is called the projected correlation function, Wp{a) 

poo 

Wp(fT) = 2 ^(cT,7r)d7r (7) 
Jo 

In our case we take the upper limit of the integration to 
be equal to Umax = 70h^^Mpc. This limit has been chosen 
to minimise the effect of the small-scale peculiar velocities 
and redshift errors (da Angela et al. 2008). If we include 
very large scales, the signal will become dominated by noise 
and; on the other hand, if we take our measurements on 
very small scales then the amplitude will be underestimated. 
Now, since Wp{a) describes the real-space clustering, the last 
equation can be written in terms of the real-space correlation 
function, ^(r), (Davis & Peebles, 1983), i.e. 

Wp{a) = 2 /"""""^ ^i^^. 

J a y(r2 -a') ^ ' 

Calculating the projected cross-correlation function Wp{a) 
will help us estimate the real-space cross-correlation func- 
tion, ^(r). The ro from the Wp{a)/a fits is estimated using 
the following equation: 

Wpja) _ /ro V r(i)r(^) 
a -UJ r(i) 

where r(.x-) is the Gamma function. 

Fig. 10 shows our Wp{a)/a measurements from the dif- 
ferent QSO samples cross-correlated with the spectroscopic 
2SLAQ LRGs. Filled circles show the results using 2SLAQ 
QS0S-2SLAQ LRGs, open circles using 2QZ QS0s-2SLAQ 
LRGs and triangles SDSS QS0s-2SLAQ LRGs. We see that 
the semi-projected cross-correlation function confirms our 
results in redshift-space, i.e. the measurements are in agree- 
ment regardless of the luminosity of the QSO sample. This 
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Figure 6. 2SLAQ QSO-2SLAQ (spectroscopic) LRG rcdsliift-spacc cross-correlation results. Filled circles show the results when the 
fibre collision effect is not taken into account and open circles show the results when the fibre collision is taken into consideration. As 
we can see the eflfect is bigger on small scales (i.e. s ^ 3h~^Mpc) than it is on larger scales. 
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Figure 7. QSO-2SLAQ LRG rcdshift-spacc cross-correlation results. Filled circles show the results when using 2SLAQ QSOs, open 
circles using 2QZ QSOs and triangles using SDSS DR5 QSOs (0.35 ^ 2: ^ 0.75, in all cases). All measurements have been made with 
spectroscopic 2SLAQ LRGs. The lines show the ^(r) fits from the photometric samples, which appear to be in agreement with the 
spectroscopic results. 
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Figure 8. §(r) fits via Limber's formula following Phillipps et al. 1977, of our w(8) measurements from spectroscopic QSO samples with 
photometric 2SLAQ and AAOmoga LRG samples. 




Figure 9. Comparison of ^(s) measurements for QSO-LRGs with QSO-QSO (da Angela et al. 2008) and LRG-LRG (Ross et al. 2007) 
measurements. Triangles are the 2QZ+2SLAQ QSO ^(s) results (0.3 < z < 2.2), filled circles show the 2SLAQ LRG ^(s) results and 
open circles show the results for the 2SLAQ QSO-2SLAQ LRG redshift-space cross-correlation. The solid line shows our fit to the 
data from 5 - 25h-iMpc, which gives so = 8.2 ± O.lh'^Mpc and 7 = l-Stov 
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Figure 10. The scmi-projcctcd cross-correlation function results for the 2SLAQ QSOs-2SLAQ (spectroscopic) LRGs (filled circles), the 
2QZ QSOs-2SLAQ (spectroscopic) LRGs (open circles) and SDSS QSOs-2SLAQ (spectroscopic) LRGs (triangles). We have also plotted 
the fits from the w{0) measurements of the photometric 2SLAQ LRG sample, using Limber's formula. 



1000 F 




Figure 11. The semi-projected correlation function results for the (2QZ+2SLAQ) QSO (open circles) and the 2SLAQ LRG-LRG 
(triangles) from da j4ngela et al. (2008) and Ross et al. (2007), respectively. The solid line shows our fit to the data from 5 — 25h~^Mpc, 
which gives ro = G.Stl'^b.-'^Mpc and 7 = l-7tl'l- 
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can also be confirmed by the fits to these measurements 
shown in Table 5. As for the case, we also include 

the fits from the w{9) measurements of the photometric 
2SLAQ LRG sample, using Limber's formula. The photo- 
metric fits are, again, in good agreement with the spectro- 
scopic measurements, further supporting the idea that the 
cross-clustering is independent of QSO luminosity. 

Fig. 11 shows the Wp{a)/a results for the 2SLAQ QSO- 
LRGs (filled circles). As in the previous Section, we have 
also included the scmi-projcctcd correlation function re- 
sults for the (2QZ-I-2SLAQ) QSO (open circles) and the 
2SLAQ LRG-LRG (triangles) from da Angela et al. (2008) 
and Ross et al. (2007), respectively. We note that, at small 
scales (^ Sh^-'Mpc), although the results are noisier than 
for the ^(s) measurements, QSO-LRG and QSO-QSO mea- 
surements have a slightly smaller amplitude than the LRG- 
LRG one but to a much lesser degree than in the ^(s) mea- 
surements. This confirms our previous interpretation that 
the amplitude difference in the redshift-space measurements 
at small scales is due to the QSO redshift errors, an ef- 
fect which does not affect the Wp{a)/a measurements. On 
larger scales, Wp(a)/a measurements confirm our previous 
observations for lower QSO-QSO amplitude comparing with 
the QSO-LRG and the LRG-LRG amplitude. Finally, the 
solid line shows our fit to the QSO-LRG data from 
5 — 25h~^Mpc (for consistency reasons with the ^(s) fits), 
which gives ro = 6.8t55:3h~^Mpc and 7 = l-7tm. This is 
similar to the (2SLAQ) LRG-LRG auto-correlation ampli- 
tude (ro = 7.45 ± 0.35h~'^Mpc) and both are higher than 
the (2QZ-I-2SLAQ) QSO-QSO amplitude (ro ~ S.Oh'^Mpc) 
at z = 1.4. 



• 2SLAQ QSO-LRG 
° 2QZ QSO-LRG 

* SDSS QSO-LRG 




Figure 12. The real-space cross-correlation function, ^(r), results 
for our different samples. The dashed lines show the fits from the 
QSO-photometric LRG w{0) measurements. 



6 THE REAL-SPACE CROSS-CORRELATION 
FUNCTION 

Using the results from the projected cross-correlation func- 
tion, described in the previous Section, and following Saun- 
ders et al. 1992, we can calculate the real-space cross- 
correlation function, £,{r), as follows: 

^W = --/ ^ da 10 

and assuming a step function for Wp{a) = Wi we finally get. 



„,,^-l^m,„r--^V'''«-'-) ,11) 

», + ^/»!-»,! 

for r = ai. 

The QSO-2SLAQ (spectroscopic) LRG real-space re- 
sults are shown in Fig. 12. In the same Figure we have also 
plotted the ^(r) fits from the QSO-photometric LRG w{6) 
measurements, described in Section 4. All the samples seem 
to give consistent results although, as already mentioned, 
these ^(r) measurements from the spectroscopic samples are 
very noisy and no significant conclusions can be drawn. Fi- 
nally, Table 6 shows the ro and 7 values from the fits to the 
spectroscopic samples, on scales of 5^ r ^ 25h~^Mpc. 




10 

a (h"' Mpc) 



Figure 13. A comparison between our 2SLAQ QSO-2SLAQ LRG 
^(<T, 7r) (solid line) results and the results using our Model I 
(dashed line). As we can see, model I is in very good agreement 
with the data both on small and large scales. 
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7 CONTRAINTS ON /3 FROM 

REDSHIFT-SPACE DISTORTIONS 

7.1 The ^(cr, 7r) cross-correlation function 

As already noted, QSO-LRG velocity dispersion and their 
dynamical infall into higher density regions are two mecha- 
nisms that distort the spherically symmetric clustering pat- 
tern in real space. We measure the ^(a, tt) cross- correlation 
function in the case of 2SLAQ, 2QZ and SDSS QSOs with 
the (spectroscopic) 2SLAQ LRGs. The purpose is to use 
the shape of the $(cr, tt) and our previous ^ measurements 
in order to model the redshift-space distortions and put 
constraints on B. Fig. 13 shows a comparison between our 
2SLAQ QSO-2SLAQ LRG ^{a, n) (solid line) results and the 
results using our model I (dashed line, see below). Model 1 
is in very good agreement with the data both on small and 
large scales. 



7.2 Description of ^{a, tt) models 

In order to place constraints on the infall parameter we 
model the redshift-space distortions measured above. First, 
we define our models for bias, 6, and infall parameter /3, 



bqbL 



(12) 



where bg and 6z, are the QSO and LRG bias, respectively, 
given by 



(13) 



where (3q and /Jl are the QSO and LRG infall parameter 
respectively and Qm is given in a flat universe as 

n^(i + 2)3 



no,(i + 2)3 + no 



(14) 



In our analysis in this Section we shall use two models. 
Model I is the one used in da Angela et al. (2008) modified 

accordingly to match our cross-correlation analysis, instead 
of the auto-correlation (QSO-QSO) which was used in that 
study. 

In general, the power spectrum in real and redshift 
space is given by (Kaiser 1987) 



Ps{k) = {l+fi{z)txlfPr{k) 



(15) 



and the similar relation between ^(r) and ^(s) is (Hamilton 
1993) 



e(r) 



(16) 



where Ps(fe) is the power-spectrum in redshift-space, Pr{k) 
is the power-spectrum in real-space, and is the cosine of 
the angle between the wavevector k and the line-of- sight. 
Equation 15 can also take the form. 



l + |/3(^) + ^/3(^)'j^o(r)Po(M) 

-^f3{z) + ip{zf^ C2(r)P2(/x) 



(17) 



H is now the cosine of the angle between r and tt, Pi{^i) 
are the Legendre polynomials of order I and Co(t'), C2(r) 
and ^4(r) being the monopole, quadrupole and hexadecapole 
components of the linear ^(r). 

In our analysis there are two infall parameters, /3q and 
Pl, one for QSOs and one for LRGs. Therefore, equations 
16, 15 and 17 (Y. P. Jing, priv. communication) should be 
modified as follows: 



1 + Ui3q{z) + l3Liz)) + II3q{z)Pl{z) 



Ps(k) = (l+/3Q(«)4)(l + /3^(«)Mfc)P.(fc) 



(18) 



(19) 



e(fT,7r) = [l + l[|3Q{z)+|3L{z)] + \|3Q{z)pL{z)^io{r)Po{^l) 
-|[/3q(z) + /3i(z)] + tjfiQ{z)[iL{z)^ 6(r)P2(/x) 

+^pQ{z)(3L{z)Wr)PA{ti) (20) 

The problem with this formalism is that the model only 
constrains the sum of the infall parameters, i.e. Pq + Pl, 
and not to each of them individually. So, in what follows we 
keep Pl constant using the value found by Ross et al. (2007) , 
Pl = 0.45±0.05 and constrain the QSO infall parameter Pq. 

The magnitude of the elongation along the 7r-direction 
of the ^{a, tt) plot caused by the peculiax velocity of the 
object is denoted by (uf)^^^, which can be expressed by a 
Gaussian (Ratcliffe et al. 1996), as 



1 



2^(a;|)i/2"^'^^ 2(c.i)i/2^ 



(21) 



To include the small scale redshift-space effects due to the 
random motions of galaxies, we convolve the ^(ct, tt) model 
with the peculiar velocity distribution, given by equation 21. 
Then, ^(ct, tt) is given by 



r 

J — c 



^'(cr, n-uj,{l + z)/H{z))f{oj,)doj. 



(22) 



where {a,TT — ujz{l + z)/H{z)) is given by equation 17. 

Finally, we exploit the Alcock— Paczynski effect ( Al- 
cock & Paczynski 1979) which says that the ratio of observed 
angular size to radial size varies with cosmology (isotropic 
clustering). If we assume that the cluster is isotropic, then 
we can constrain the cosmological parameters by requiring 
that they produce equal tangential and radial sizes. In our 
case the angular and the radial size are a and n. Keeping 
the meanings of test and assumed cosmology as used by da 
Angela et al. (2008) and following their fitting procedure we 
obtain our results (see below). In particular, we use the 7 
values from our fits to ^(s), let ro and the velocity dispersion 
vary as free parameters and compute the values for each 

- p{z) pair. 

The second model (model II) is as follows, ^{a, tt) is 
now deflned as (Peebles 1980, Hoyle 2000) 



1 + 



a, tt)= f 

J —c 



{1+ ar))f {oJ.)doj. 



(23) 



where the /(wz) is given, as before, by equation 21. Next 
we introduce the infall velocity of the galaxies, v{rz), as a 
function of the real-space separation along the tt direction. 



QSO-LRG 2-Point Cross-Correlation Function and Redshift- Space Distortions 13 



Vz- This can be derived from the equation of conservation 
of particle pairs, which are within a comoving separation r 
from a mass particle, i.e. (Peebles 1980) 

^ rx'r(^c,t)dx+ J^r^(l+r(r,t))w(r,t) = (24) 
Jo 



St 



a(t) 



where a{t) is the scale factor. Assuming that ^(r) is de- 
scribed by a power law model, we solve the above equation 
to find the infall velocity of the particles 



vir.) = -^fi.(.)»-«i/(.)r.— if^^ 
3-7 bQbL+^QL{r) 



(25) 



Wo now modify equation 23 to include the effects of the bulk 

motions 

/ + 00 
{l + i{r))f{uj4l + z)-v{rz))dcj. (26) 
-oo 

This is the model II ^(a, tt); we then follow the same 
implementation of the "Alcock-Paczynski" effect and fitting 
procedure as for model I. As in model I, we keep 6i constant, 
using the same value as above (Ross et al. 2007), i.e. /?i = 
0.45 ± 0.05 and 6l = 1.66 ± 0.35. 



7.3 Results 

We now use our ^{a,n) measurements from the previous 
Section and the so and 7 values from the ^(s) fits shown 
in Table 4 to put constraints on /3q and bq and {uj^)^^^ for 
each of the QSO samples. The results are shown in Table 7 
and in Figures 14-19. 

Comparing the results from both models from the three 
different QSO samples associated with the same spectro- 
scopic 2SLAQ LRG sample, we note that model I gives 
slightly lower QSO-LRG velocity dispersions than model II, 
but both are consistent (~620kms~^ and ~ 727kms~^) with 
the expected (wj)^''^ ~728kms^^, produced by adding in 
quadrupole QSO and LRG velocity dispersions from previ- 
ous QSO-QSO and LRG-LRG studies, which gave SOOkms"^ 
and 330kms^^, respectively. Comparison of bq between the 
different samples (for both models) shows that the results 
are in good statistical agreement (Table 7). The best way 
to summarise the results is to average them. All six mea- 
surements (the three samples and the two models) give an 
average of /3q = 0.55 ± 0.10 and bq = 1.4 ± 0.2 at z=0.55, 
which is consistent with the values found by da Angela et 
al. (2008), Pq = O.eotg.n and bq = 1.5 ± 0.2 at z=1.4 and 
bq = 1.1 ± 0.2 at 2; ~ 0.6 (see their Fig. 13). 



8 QSO BIAS AND HALO MASSES 

8.1 QSO-LRG clustering dependence on 
luminosity 

Previous attempts to study the dependence of clustering 
on luminosity, were not successful because of the redshift- 
luminosity degeneracy, as higher luminosity QSOs reside 
at higher redshifts. Here, combining QSOs with large LRG 
samples we get the statistical power to break this degener- 
acy. In this Section, we shall try to examine if and how the 
QSO-LRG clustering depends on QSO luminosity, at fixed 
redshift. 



We first estimate the average absolute magnitude of 
each of our QSO samples, as follows: 

Mbj = 6j - 25 - 5Iogiodt -|- 2.5{1 + a')logio(l + z) (27) 

where Mbj is the absolute magnitude of each QSO, bj (or 
g) is its apparent magnitude , di, is the luminosity distance 
(Mpc) that corresponds to the redshift z and the last term 
is the k-correction, where we have assumed a QSO spectral 
index a' = —1.0. We should note here that we treat the 
bj band (2QZ) as being equivalent to the g band (2SLAQ, 
SDSS; Richards et al. 2005). 

To check if there is a dependence of QSO-LRG cluster- 
ing on luminosity wc use the integrated correlation function, 
as a more robust statistical tool (see da Angela et al. 2008). 
We calculate it up to scales of 20h~^Mpc and normalise the 
result to the volume contained in a sphere with radius of 
20h-^Mpc: 



C20 = 



203 



f''" 2 
/ C(s)s ds 

Jo 



(28) 



The choice of the radius has been made for two reasons. 
The 20h~^Mpc scale is large enough to apply linear theory 

(Croom et al. 2005), and redshift-space distortions (finger 
of god or redshift errors) do not significantly affect the mea- 
surements. 

In the case of the QS0s-2SLAQ (spectroscopic) LRGs 
we have estimated ^20 via equation 28 using the ^(s) mea- 
surements shown in Fig. 7. For the QS0-2SLAQ (photo- 
metric) and AAOmega LRGs we have substituted the ^(s) 
in equation 28 with the ^(r) fits shown in Fig. 8. The re- 
sults using 2SLAQ LRGs (spectroscopic and photometric) 
are shown in Fig. 20 and using AAOmega LRGs are shown 
in Fig. 21. No conclusion can be drawn about the redshift 
evolution of the QSO-LRG clustering as the average red- 
shift of the samples is too restricted. Although the ^20 re- 
sults using the AAOmega LRGs, show some indications that 
bright QSOs (SDSS) cluster less with LRGs than faint QSOs 
(2SLAQ), the results using the 2SLAQ LRG samples (both 
spectroscopic and photometric) stay statistically constant, 
thus confirming the results in the previous Sections, that 
the QSO-LRG clustering is independent of QSO luminosity. 

8.2 QSO bias 

Having tested the luminosity dependence of the QSO-LRG 
clustering, we now investigate the dependence of the QSO 
bias on the luminosity. Assuming that this bias is indepen- 
dent of scale, we calculate the QSO bias using: 



bqbi 



1 ^QLjr, 20) 
bL^ mm (r,20) 



(29) 



where is the matter real-space correlation function, aver- 
aged in 20h~^Mpc spheres. In the case of QSO-photometric 
LRGs we use our ^(r, 20) measurements as estimated before. 
For consistency, for the QSO-2SLAQ (spectroscopic) LRGs, 
wc use our Wp{a) measurements (shown in Table 8) since 
they are less noisy than those for ^(r). To estimate ^mm we 
use the values as estimated by da Angela et al. 2008. So in 
the case of QSO- AAOmega LRGs ^„„„ = 0.11 at {z ^ 0.7), 
and in the case of QSO-2SLAQ LRGs ^mm = 0.12 {z « 
0.55). Finally, we calculate the LRG bias for the AAOmega 
and 2SLAQ, based on the results shown in Table 4 of Ross et 
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al. (2007) (AAOmega: ro = 9.03 ±0.93 and 7 = 1.73 ±0.08, 
2SLAQ: ro = 7.45 ± 0.35 and 7 = 1.72 ± 0.06). We find, 
bL(AAn) = 2.35 ±0.20 and bL(2SLAQ) = 1.90 ±0.08. The lat- 
ter is in reeisonable agreement with the value found by Ross 
et al. (2007), from rodshift-spacc distortions, 6_l(25laq) = 
1.66 ± 0.35. The derived QSO bias values for each case are 
shown in Table 8 (as well as the corresponding /3q values) 
and in Figures 22 and 23. 

Comparing the values for the QSO biases from the dif- 
ferent samples, wo note that the QSO biases using 2SLAQ 
LRG samples show indications for luminosity dependent 
QSO bias, in the sense that bq reduces for higher luminosity 
samples, at least in the case of spectroscopic 2SLAQ LRGs. 
The same pattern is repeated when using AAOmega LRG 
samples. The spectroscopic samples yield lower bq values 
than the photometric samples. This is due to the fact that 
the amplitude of the ^(r) measurements of the photometric 
samples is higher than the amplitude of the it)p(cr) measure- 
ments of the spectroscopic samples (see Fig. 10). Combining 
the 2SLAQ (photometric and spectroscopic) samples with 
the photometric AAOmega samples we find 6q = 1.90±0.16, 
bq = 1.85 ±0.23 and 6q = 1.45 ±0.11, for 2SLAQ, 2QZ and 
SDSS QSOs, respectively. Comparing now the values for the 
QSO bias from the spectroscopic 2SLAQ LRG samples with 
those obtained in Section 7.3, we note that the amplitude re- 
sults, give an average of 6q = 1.5 ±0.1 which is in very good 
agreement with the average of bq — 1.4 ±0.2 obtained from 
the redshift-spacc distortion results. Our measurements give 
an overall QSO bias of bq « 1.5 aX z = 0.55 and Mbj k, —23. 

In Figures 22 and 23 we have also plotted two points 
(stars), that are at low redshifts, taken from Fig. 13 of da 
Angela et al. (2008). The one with Mij ~ -24.0 at a ~ 
0.7 is in statistical agreement with our bq values from the 
AAOmega LRG samples, at the same mean redshift and 
brightness. The second one with M^j ~ —22.9 at ^ ~ 0.6 
is lower than our hq values from 2SLAQ LRG samples, at 
z = 0.55, but statistically not rejected by them (at least 
not by those from the spectroscopic samples). The overall 
impression is that our bq « 1.5 aX z = 0.55 is in agreement 
with the values found by da Angela et al. (2008), bq = 
1.5 ± 0.2 aX z = 1.4 and slightly higher than bq = 1.1 ± 0.2 
found aX z ^ 0.6. 



8.3 Dark Matter Halo Mass 

Since the bias of Dark Matter Halos is correlated to their 
mass (Mo & White 1996), we shall attempt to measure this 
mass (Mdmh)- In our analysis we shall follow da Angela 
et al. and Groom et al. and assume an ellipsoidal collapse 
model, described by Sheth et al. (2001). The bias and the 
Mdmh are related via 



is the rms fluctuation of the density flcld on the mass scale 
with value Mdmh and G{z) is the linear growth factor (Pee- 
bles 1984). The a{MDMH) can then calculated as 



b{ML 



f) = 1 + 



1 r 0.5/ 2., . 0.5./ 2.,1-c 

[q (af ) + a 0(01/ ) 



aO-^Sc{z) 

(av-y 



ai^^ + b{l-c){l-^y 



(30) 



where a — 0.707, b — 0.5 and c — 0.6. f is defined as 
f = 5c{z) / (j{Mdmh , z), with 5c to be the critical density for 
collapse, given by, 5^ = 0.15(127r)ir2„(2;)° °°^^ (Navarro et 
al. 1997). a(MDMH,z) = a(MDMH)G(z), where a{MDMH) 



(31) 



1 

<7{MDMHf = ^ j k'P{k)w{krfdk 

with P(k) to be the power spectrum of density perturbations 
and w{kr) is the Fourier transform of a spherical top hat, 
which is given by (Peebles 1980): 

, , sin(A;r) - A;rcos(A;r) , , 

^'(fcO = 3 ^3 (32) 

where the radius and mass are related through 

where po is the present mean density of the Universe, given 
by Po = f^mPcHt = 2.78 X 10^^Ql,h'^MQMpc-^. The power 
spectrum used in our analysis has the linear form, P{k) = 
PoT{k)^kP' , with Po to be a normalisation parameter which 
depends on as and T(k) is the transfer function (Bardeen et 
al. 1986). 

The results arc shown in Figures 24 and 25. Once again, 
although for the AAOmega LRG samples the derived QSO 
halo masses show indications of increasing as we move to 
fainter QSO samples, in the case of 2SLAQ (photometric and 
spectroscopic) LRG samples halo masses stay statistically 
constant. The average value is Mdmh = IQ^'^hT^ Mq. Com- 
paring now this result with those from other authors (Groom 
et al. 2005, da Angela et al. 2008), we note that their Mdmh 
estimates are generally lower than ours (~ 3 x IQ^^h^^ Mq) 
although at higher redshifts {z = 1.4). They also find that 
the hosts of QSOs have the same mass at all redshifts, thus 
rejecting cosmologically long-lived QSO models. Our higher 
masses aX z = 0.55 may be more consistent with the long- 
lived predictions of 6x 10"/r^MQ at z = and 10"/r 
at 2; ~ 0.5. The caveat is that for our measurements we need 
to use a value for bh in order to derive bq . 



9 DISCUSSION + CONCLUSIONS 

In this paper we have performed an analysis of the clustering 
of QSOs with LRGs. For this purpose, we first used the 2- 
point angular cross-correlation function, w{6), and measured 
the cross-correlation between 2SLAQ and AAOmega LRGs 
and different luminosity QSOs. The results show that there 
is little cross-correlation dependence on QSO luminosity. 

Next, we measured the redshift-spacc cross- correlation 
function. We again see no QSO-LRG clustering dependence 
on QSO luminosity, as all the QSO-spectroscopic LRG sam- 
ples gave similar results. We used Limber's formula to fit ro 
to 2-D results. The fits for ro from 3-D ^(s) are in very good 
agreement with the fits to the 2-D w(6) results. Then, we 
compared our QSO-LRG clustering with 2SLAQ LRG-LRG 
(Ross et al. 2007) and 2QZ+2SLAQ QSO-LRG (da Angela 
et al. 2008) clustering results. On small scales, the QSO-QSO 
and QSO-LRG results appear flatter than the LRG-LRG re- 
sults. As confirmed later by the Wp{(j)/a measurements, this 
appears to be due to the larger QSO redshift errors (broad- 
lines). On larger scales (^ 5h~^Mpc) the QSO-QSO cor- 
relation amplitude appears lower than the QSO-LRG and 
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LRG-LRG amplitudes, suggesting that 6q < Bl- The frac- 
tional errors on ^ql are ~ 50% smaller than those on ^qq 
in same redshift range. 

The results from the semi-projected cross-correlation 
function, Wp{a)/a, arc in agreement with our ^(s) obser- 
vations, yielding consistent ro and 7 values. The compar- 
ison with the 2SLAQ LRG-LRG (Ross et al. 2007) and 
2QZ-I-2SLAQ QSO-LRG (da Angela et al. 2008) confirms 
our results from redshift-space, i.e. that the small-scale am- 
plitude difference in ^(s) is due to the larger QSO redshift 
errors and that QSO-QSO clustering amplitude is lower than 
QSO-LRG and LRG-LRG amplitude. The real-space cross- 
correlation function, 5(r), also seems to agree with the ^(s) 
and Wp{(j)/a results, although it is noisier. 

Then, we measured the ^(cr, tt) cross-correlation func- 
tion for all our (spectroscopic) samples and used the results 
to model the redshift-space distortions. For that, we used 
two models which gave consistent results with each other 
and between the different samples. The redshift-space dis- 
tortions yielded an average of Pq = 0.55±0.10, 6q — 1.4±0.2 
which is slightly higher than 6q = 1.1 ±0.2 at a ~ 0.6, from 
da vlngela et al. (2008). Note that this latter result does not 
come from analysing the amplitude of ^qq; there being too 
few QSO pairs to make redshift-distortion auto-correlation 
analysis viable at z ~ 0.6. 

After calculating the average absolute magnitude of 
each QSO sample we measured the integrated cross- 
correlation function for each one of our QSO-LRG samples. 
No evidence was found for QSO-LRG clustering dependence 
on QSO luminosity. Then, using LRG biases as estimated 
by previous studies, we calculated the QSO biases. There 
were indications of luminosity dependence, in the sense that 
bq may reduce as we move to brighter QSO samples. Our 
analysis yielded a. bq = 1.5 ± 0.1 (at Mi,j w —23) which is 
in very good agreement with our result from redshift-space 
distortions. 

Finally, using the relation between bias and Mdmh sug- 
gested by Sheth et al. 2001, we calculated the correspond- 
ing masses of the QSO hosts. QSO halo masses were esti- 
mated to be ~ 10" /i~^Mq at « « 0.55. Our Mdmh esti- 
mations arc higher than those from Groom et al. (2005) and 
da Angela et al. (2008) (at z — lA) and may be explained 
by long-lived QSO population models. Since the bias val- 
ues are independent of QSO luminosity at fixed redshift, 
the halo masses are also independent of luminosity and this 
represents the main result of this paper. 
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Figure 14. Likelihood contours of - l3{z = 0.55) for 2SLAQ 
QSO— 2SLAQ LRGs, using model I. A ACDM cosmology is as- 
sumed along with a model where ro = 8.2, 7 = 1.6 with (w^)^^"^ = 
630kms-^ The best fit value is /3q{z = 0.55) = 0.45to;22. 
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Figure 17. Likelihood contours of Cl^ - P{z = 0.55) for 2SLAQ 
QSO— 2SLAQ LRGs, using model IL A ACDM cosmology is as- 
sumed along with a model where ro = 8.2, 7 = 1.6 with (w^)^^^ = 
750kms-i. The best fit value is /3q(2 = 0.55) = 0.40lo;22. 
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Figure 15. Likelihood contours of — f3{z = 0.55) for 2QZ 
QSO-2SLAQ LRGs, using model I. A ACDM cosmology is as- 
sumed along with a model where ro = 8.0, 7 = 1.7 with {w'j.)^^^ = 
560kms-i. The best fit value is I3q{z = 0.55) = 0.55lo:||. 
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Figure 18. Likelihood contours of - l3{z = 0.55) for 2QZ 
QSO-2SLAQ LRGs, using model II. A ACDM cosmology is as- 
sumed along with a model where ro = 8.0, 7 = 1.7 with (w^)^^^ = 
710kms-i. The best fit value is I3q{z = 0.55) = Omto^l- 



0.8 - 




0.2 0.4 0.6 0.8 1.0 



Figure 16. Likelihood contours of n'^ - I3{z = 0.55) for SDSS 
QSO-2SLAQ LRGs, using model I. A ACDM cosmology is as- 
sumed along with a model where ro = 7.5, 7 = 1.8 with (w^)^^^ = 
670kms-i. The best fit value is /Sgiz = 0.55) = O.GOto.H- 
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Figure 19. Likelihood contours of - /3{z = 0.55) for SDSS 
QSO-2SLAQ LRGs, using model II. A ACDM cosmology is as- 
sumed along with a model where ro = 7.5, 7 = 1.8 with (ui^)-'/^ = 
710kms-i. The best fit value is /3(z = 0.55) = Om^ojiT- 
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Figure 20. ^20 cross-correlation measurements of the three QSO samples with 2SLAQ LRGs. Filled symbols show the results using 
spectroscopic samples and open symbols using photometric (LRG) samples. 
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Figure 21. {20 cross-correlation measurements from the three QSO samples with (photometric) AAOmega LRGs. 
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Figure 22. Measurement of the QSO bias, bq, for QSOs and 2SLAQ LRG samples. For consistency, spectroscopic samples use ^20 from 
Wp{a) in Table 8 rather than the ^(s) values shown in Fig. 20. Stars show the two points taken from Fig. 13 of da Angela et al. (2008). 
The fainter one is at (z) ~ 0.6 and the brighter at (z) ~ 0.7. 
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Figure 23. Measurement of the QSO bias, 6q, for QSOs and the (photometric) AAOmega LRG sample. Stars show the two points 
taken from Fig. 13 of da Angela et al. (2008). The fainter one is at (z) ~ 0.6 and the brighter at (z) ~ 0.7. 
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Figure 24. Measurement of the Mj^mh, for different QSO and 2SLAQ LRG samples. 
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Figure 25. Measurement of the Mumh, for different QSO and AAOmega LRG samples. 
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Table 3. ro and 7 values from the fits on the QSO-2SLAQ (photometric) LRG ^(r) measurements. 



photometric 2SLAQ LRGs 



2SLAQ QSOs 2QZ QSOs SDSS QSOs 

'^O '-^-O.S °-"-0.4 '-"-0.3 

1 -+0.2 , -+0.2 , jj+0.1 



Table 4. so and 7 values from the fits on the QSO-2SLAQ (spectroscopic) LRG ^(s) measurements, on scales of 5-25h ^Mpc. 



spectroscopic 2SLAQ LRGs 



2SLAQ QSOs 2QZ QSOs SDSS QSOs 



so 



,+0.1 o ri+0-3 7 c+O.S 



,+0.2 , -+0.2 , 0+0.2 

■■-0.1 -'■•'-O.l -'■•°-0.3 



Table 5. ro and 7 values from the fits on the Wp{a)/<7 measurements, on scales of 5-25h ^Mpc. 



Wp{a)/a 



spectroscopic 2SLAQ LRGs 



2SLAQ QSOs 2QZ QSOs SDSS QSOs 



ro 



D + O.l 



, 7+0.2 , c+O.l -I 0+0.1 

^•'-0.3 ^•'-'-0.2 ^-"-O.l 



Table 6. ro and 7 values from the fits on the ^(r) measurements, on scales of 5-25h ^Mpc. 



spectroscopic 2SLAQ LRGs 



2SLAQ QSOs 2QZ QSOs SDSS QSOs 

ro '-"-0.1 '-"-0.3 °-'^-0.2 



,+0.2 , f.+0.1 9 0+0.2 

■^-0.1 -"-"-O.l ^•'^-0.3 
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Table 7. QSO bg, /9q and (a;^)^/^ measurements from modelling the redshift-space distortions. 



QSO 6q and Pq 




2SLAQ QSOs 


model I 
2QZ QSOs 


SDSS QSOs 


2SLAQ QSOs 


model II 
2QZ QSOs 


SDSS QSOs 




l.DD_Q gg 


1 S6+^*l 

l.CSD_Q gy 


-, r,c+0.88 

J^-^J_0.46 


1 0*7+2.28 
-'-■°'-0.83 


1 25+l-'^* 


-'^■^'^-0.32 




<-'-<iO_Q 22 


"i^+o-ss 

U.OO_Q 3g 




n 40+0-32 


0.60t°Ji 


0.65t°Jf 


(a;2>i/2 (kms-i) 


630 


560 


670 


750 


720 


710 



Table 8. QSO-LRG ^20 cross-correlation measurements, as well as QSO feg and /9q measurements, assuming 6_L(^An) = 2.35 ±0.20 and 
bL(2SLAQ) = 1-90 ± 0.08, from the amplitude results. For consistency, the ^20 measurements for the spectroscopic samples come from 
u>p((t) and from ^(r) via w(0) for the photometric cases. 



QSO bq and /3q 




spectroscopic 2SLAQ LRGs 


photometric 2SLAQ LRGs 


photometric AAOmoga LRGs 




2SLAQ QSOs 


2QZ QSOs 


SDSS QSOs 


2SLAQ QSOs 


2QZ QSOs 


SDSS QSOs 


2SLAQ QSOs 


2QZ QSOs 


SDSS QSOs 


^20 


0.37 ± 0.06 


0.33 ± 0.10 


0.31 ±0.14 


0.44 ±0.08 


0.49 ± 0.06 


0.38 ±0.08 


0.55 ±0.11 


0.50 ±0.11 


0.33 ± 0.08 


bQ 




1 44+0-26 


1 07+0.37 

-0.37 


1.911^2! 


9 1^+0.16 


1 67+0-1* 

-0.18 


9 17+0.28 

^■-^ ' -0.28 


1 95+0.27 


1 ,(,+0.19 


00.6 


46+° °^ 
"■^"-0.05 


^1+0-09 


54+0-14 


0.28iro| 


S5+0-02 
U.OO_Q 02 


45+0-OS 


PI ,7+0.05 


41+0-06 

"•^^-0.06 


n CI +0.10 

0-61Io.io 




-21.7 


-22.9 


-23.7 


-21.7 


-22.9 


-23.7 


-22.2 


-23.5 


-24.5 



